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C^*) • We show that additional p-wave component is generated in the pure d-wave superconductor in 

the presence of the uniform current. When the current flows in the antinodal direction the spectrum 
has a gap over all Fermi surface. If the current flows in the nodal direction, the gap opens only in the 
direction parallel to the current. The correction to the current due to the p-wave component of the 
order parameter is linear in vector potential. We show, that the spin-orbit coupling is responsible 
for the additional component of the order parameter. 
Recently it was proposed that classification of the order parameter in the ferromagnetic superconductor may be 
performed in accordance with corepresentations of nonunitary magnetic group [1—3]. Similar approach was used 
in Ref. [4] to consider the effect of broken time-reversal symmetry at T* on the superconducting state in high-T c 
superconductors. Here we would like to extend this type of analysis for the case of a d-wave superconductor with 
uniform current or/and the external magnetic field. 

There is a common belief that the high temperature superconductors have nontrivial order parameter transforming 
as Big (x 2 — y 2 ) representation of point group. When external magnetic field is applied the symmetry of the 
O order parameter is reduced [5,6]. As a result, for magnetic field directed along the c-axis the new order parameter has 
the form: d x i_ y i + id xy . Furthermore, additional component of the order parameter may be induced by the surface. 
Since any surface breaks inversion symmetry, additional p wave component is generated near the surface due to the 
surface induced spin orbit coupling [7]. In the following we show, that a constant superconducting current reduces 
the symmetry of the order parameter to d x 2_ y 2 + ip x when the current is directed along the x-axis. Such form of 
the order parameter is consistent with the group theory analysis and may be generated by the spin-orbit coupling. 
We will construct the corresponding invariant terms in the free energy and show how p-wave is generated from the 
Ch j microscopic BCS theory. 

Since a constant superfluid current is a real vector that changes sign under time reversal transformation, the total 
group will be reduced to nonunitary (magnetic) group D 2 h(C2 V )- We consider for simplicity that superfluid 
current is directed along the £-axis. Unitary subgroup Civ has four elements E,C2 X ,o~h,o- y and four representations 
Ai, A2, B\, £?2- In accordance with the general theory of corepresentations [8] all 4 irreducible representations belong 
• to the class 'a'. This means that each irreducible representation generates one nonequivalent corepresentation of the 
i—\ ' nonunitary group. All corepresentations are one dimensional and are listed in the Table 1 together with the basis 
CN . functions. 
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TABLE I. Corepresentations of the nonunitrary group, when current flows in the x direction. 
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Here a and b are real numbers and e = exp (2i<f>) where <f> is real. Following Ref. [2] we introduced factor 2 in 
the exponent. To clarify the derivation of basis functions in Table 1 we construct the basis function for the A\ 
corepresentation. It is easy to see that exp (—i(f){k x — k y ) transforms as the A\ representation of the unitary subgroup 
as well as corresponding corepresentation of the nonunitary group. Since the group in the presence of uniform current 
' does not have the center of inversion, a mixture of singlet and triplet order parameter is possible. The general 
expression for the triplet part of the order parameter has the following form: 

*i = xf x (k x ,k y ) + yfy{k x ,k y )-\-zf z (k x ,k y ) (1) 

where x, y, z are spin unit vectors in the a, b, c axis directions respectively. The functions f x . VtZ {k Xl k y ) are odd with 
respect to the k — > — k transformation. Applying the ah operation to f x ^ y (k x , k y ) we obtain that f x , y (k x , k y ) = 
—f x .y(k x , k y ) for A\ and B\ representations. Therefore f XlV (k x , k y ) = for A\ and B\ cases. Applying all operations 
from Table 1 to the remaining function f z (k x , k y ) we obtain the following set of equations: 

fz{k x ,—k y ) = — f z (k x , k y ) 



1 



fz{ k x , ky) — ef z (k Xl k y ) 
fz{~k X i ky) = —ef z (k Xl ky) 



(2) 



It is easy to see that the function 

®(k x , k y ) = exp {-i4>)(d x 2_ y2 + it]p x ) (3) 

where p x — zk y satisfies all the equations and transforms as the A\ corepresentation of the nonunitary group D 2 h{C2v)- 
This means that the uniform current in the x direction generates the p x wave contribution to the gap function. The 
gap function in that case is defined as A(k x ,k y ) = i(d x 2_ y 2 + ir\uj ) ky)u 2 . Since AA^ oc <to the corresponding phase 
is unitary and the gap in the excitation spectrum is determined by \d x 2_ y 2\ 2 + rfk 2 [9]. As a result, the gap is fully 
developed over all Fermi surface. 

Let us show that the term in the free energy: 

iV ( v xPx - v y p y )d* x 2_ y 2 + c.c, (4) 

is an invariant of the group D±h- Here r\ is a real constant and (v x ,v v ) are components of the superfluid velocity. 
Indeed, (v x ,v y ) and (p x ,p y ) transform as the E u representations of the group D^. Since the direct product E u x E u = 
Ai g + A 2g + Big + B 2g contains B lg (d x 2_ y 2) representation, Eq.(4) represents the true scalar. The presence of i in 
Eq(4) is important since the superfluid current is antisymmetric with respect to the time reversal symmetry. 

The existence of the p component of the order parameter in the presence of the supercurrent may be seen if we 
write relevant Lifshitz invariant in the free energy [10]: 

F = F d + iri\p d {D xPx - D vPy ) - MKPl ~ D* y p* y )) + a p (\ Px \ 2 + \ Py \ 2 ) (5) 

where D\ = — iVj — 2eA; (I = x, y, z), A is a vector potential, and Ed is the free energy for the <i-wave superconductor 
without curren. Assuming that tpd = ipdoexp{imvx) and (p x ,p y ) = (p x ,0)exp(imvx) we obtain a contribution to 
the free energy, which is similar to Eq.(4). By minimizing Eq.(5) with respect to p x we obtain the amplitude of the 
p-component of the order parameter: p x — ir) mvipd/otp. 

When the current is flowing in the nodal direction, Bi g (x 2 — y 2 ) representation of the D±h group is no longer A\ 
corepresentation of the reduced nonunitary group. Group theoretical analysis for this case is presented in Table 2. It 
should be pointed out that this case is similar to the case II considered in the Rcf. [4] . 

TABLE II. Corepresentations of the nonunitrary group, when current flows in the nodal direction. 
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Similarly to the case when the current flows in the x direction, the corresponding phase is unitary: AA^ oc ao, and 
the gap function has the form: \d x 2_ y 2\ 2 + a 2 (k x + k y ) 2 . Contrary to the previous case, however, gap opens only in 
the direction parallel to the current, and the spectrum remains gapless in the direction perpendicular to the current 
[11]. 

Let us now discuss the correction to the current when the additional component of the order parameter is generated. 
It is easy to see from Eq.(5) that the additional contribution to the current is given by: 

jx = jdx + 2ieri (ip d Px ~ tpdP* x ) (6) 

Substituting p x = ir] mvipd/a p to the last equation and taking into account that in the quasiclassical limit mv = —2eA, 
correction to the supercurrent has the form: 

j =id + ^ 2 7 1 " 2 A\ip d \ 2 /a p (7) 

At this point we should point out that usually the nonlinear London equation is discussed in the case of the d-w&ve 
superconductors with the nodes in the spectrum [11]. Nonlinear corrections to the London equation appear due to the 
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Dopplcr shift of the spectrum in the presence of the superfiow. As a result, in some regions of the Fermi surface the 
excitation energy becomes negative leading to the finite quasiparticle current. However, as we have shown here, the 
generation of the p-wave component of the order parameter in the presence of the superfiow leads to the correction 
to the supercurrent which is linear in the vector potential. Opening of the p wave component of the order parameter 
leads to the linear in current correction to the penetration depth and could be detected experimentally. 

At the end we suggest one possible microscopic mechanism, which causes the appearance of the p-componcnt of 
the order parameter in the presence of uniform current. We assume that current flows along the x axis. Spin-orbit 
coupling in that case could be written as 

H so = !7^^[v x 9 k ]^cT3^k 

k 

where v is the superfluid velocity and 7 is the spin-orbit coupling constant. Following Balatsky [5] we can calculate 
the correction to the anomalous Green function F(k, u) = A (k)/D(k, u), where D(k,u) = u 2 + £(k) 2 + |A (k)| 2 , 
using perturbation theory due to the interaction H so . The first correction to F°(k, u) is: 



SF(k,u) = z 7 «G (k, W )d fey a 3 F°(k,^) 
(i u - g(k)) c7 3 n . 
D(k,u) 2 
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(— iu + £(k))«Ti . . 

= D(ku)" sm[ky) (8) 

Here we take into account that A (k) = A (cos(k x ) — cos(k y ))ia2, and G°(k,u) = (iu — £(k))<ro/.D(k, u). To 
estimate the p-wave correction to the order parameter we assume a repulsive separable interaction in the following 
form: V y (k, k ) = V v sin (k y ) sin (fc ) As a result, the correction to A (k) is given by: 

Ai(k) =TJ2V V sin (ky) sin (k' y )SF(k , u) 

A 

= -iC-yv—NoVy sin (k y )a x (9) 

Here C is a real constant of the order of 1, and N is the density of states at the Fermi surface. 

In conclusion, we have shown, that additional p wave component appears in the case of d-wave superconductor 
in the presence of uniform supercurrent. This effect leads to the additional contribution to the current linear in the 
vector potential. It is shown that spin orbit coupling is responsible for this effect. 

I am grateful to A.S. Alexandrov, A. Balatsky, D.F. Agterberg, D. Mihailovic, T. Mertelj, and J. Demsar for useful 
discussions. 



[1] LA. Fomin, Pis'ma ZhETF, 74, 116 (2001). 

[2] LA. Fomin, ZhETF, 122, 1089 (2002). 

[3] V.P. Mineev, cond-mat/0204263. 

[4] R.P Kaur and D.F. Agterberg, cond-mat/0306677. 

[5] A.V. Balatsky, Phys. Rev. B,61, 6940 (2000). 

[6] T. Koyama and M.Tachiki, Phys. Rev. B,53, 2662 (1996). 

[7] V.M. Edelstein, Phys. Rev. Lett., 75, 2004 (1975; L.P. Gorkov, E.I. Rashba, Phys. Rev. Lett., 87, 037004 (2001). 
[8] V.V. Eremenko, 'Vvedenie v opticheskuiu spektroskopiu magnetikov,Naukova dumka, Kiev, (1975). 
[9] M. Sigrist and K. Ueda, Rev. Mod. Phys., 63, 239 (1991). 
[10] V.P. Mineev and K.V. Samokhin, JETP, 78, 401, (1994). 

[11] We have to pint out that the quasiparticle spectrum has additional term due to the Doppler shift. As a result spectrum 
could be negative in some regions in the Fermi surface. V.P. Mineev, K.V. Samokhin, Introduction to unconventional 
superconductivity.; Gordon and Beach publishers, 1999. 



3 



